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We study superconductivity in Boron using first-principles LAPW calculations, the rigid-muffin- 
tin approximation and the McMillan theory. Our results point to an electron-phonon mechanism 
producing transition temperatures near lOK at high pressures in agreement with recent measure- 
ments. 
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I. INTRODUCTION 

In a recent paper Eremets et a/.lll reported experiments 
where boron transforms from a non-metal at normal pres- 
sures to a superconductor at very high pressures above 
160 GPa. They presented results that showed an in- 
crease of the superconducting temperature, Tc, from 6K 
at 175 GPa to IIK at 250 GPa. 

Previous theoretical studies on the metallization of 
boron based on first-principles total-energy calculationscl 
predicted that the 12 atom insulating form of boron 
(q;i2B) at atmospheric pressure transformed to a metallic 
body-centered tetragonal phase at 210 GPa and subse- 
quently to an fee structure at 360 GPa. 

In this work we performed total energy and band struc- 
ture calculations for the rhombohedral 0:126 and fee 
phases. Our results for the fee phase were used as an in- 
put to the rigid-muffin-tin (RMT) theory of Gaspari and 
GyorffyO to determine the value of the Hopfield param- 
eter T] and subsequently evaluate the coupling constant 
A and Tc. Our results for the parameter give values 
comparable to those reported in the pasta for metallic 
hydrogen at very high pressure. Assuming that the fee 
phase is the real high pressure structure in the experi- 
ments of Eremets et alB, this indicates that the simple 
RMT theory accounts for high-pressure superconductiv- 
ity in boron. 



II. THEORY AND RESULTS 

In this work we have applied the RMT theoryf3 to cal- 
culate the Hopfield parameter 77 given by the expression: 
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where Si is the scattering phase shift at the Fermi energy 
Ep and angular momentum £, ivj^"* is the single-scatterer 
density of states which, as defined in Ref. ^, is an integral 
involving the radial wave functions. N{Ef) is the total 
density of states (DOS) at Ep and Ng are the angular 
momentum components of the DOS inside the muffin-tin 
spheres. Eq. (|l|) is exact to ^ = 1, but for higher values 
of £ it involves non-spherical corrections. The necessary 
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FIG. 1: The electronic band structure of ai2B, calculated by 
the LAPW method. The solid vertical lines represent high- 
symmetry points, while the dashed vertical lines represent 
zone boundaries. The horizontal line near 0.85 Ry is the Fermi 
level. 



input to Eq. (Q) was generated from a set of full potential 
Linearized Augmented Plane Wave (LAPW) calculations 
that we performed for fcc-B using touching muffin-tin 
spheres, and the tetrahedron method for the DOS. 

We also performed an LAPW calculation for 0:126. This 
is a rhombohedral structure, space greup RSm-Dg^^ 
(#166), with measured lattice parameteril a = b = c = 
9.56 a.u. and angles a = (3 — j = 58.06°. The en- 
ergy bands of 0126 are shown in Fig. |l]. We obtain 
a valence bandwidth of 15.9 eV and an energy gap of 
1.47 eV separating valence and conduction bands. The 
angular momentum character of the states can be seen 
from Fig. ^ which shows the expected dominance of the 
p-states through the whole spectrum. 

For fee B we calculated an equilibrium lattice param- 
eter a=5.37 a.u. and a bulk modulus of 282 GPa, to be 
compared with the values of 5.34 a.u. and 269 GPa re- 
ported by Mailhiot et aln We note that for a — 4.60 a.u., 
which, according to our total energy calculations corre- 
sponds to a pressure of 307 GPa, the nearest neighbor 
B-B distance is 3.25 a.u., which is close to the B-B dis- 
tance (3.37 a.u.) in the newly discovered superconductor 
MgB2. This leads us to believe that superconductivity 
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FIG. 2: The electronic density of states of 0126, computed 
by smearing the eigenvalues of an LAPW calculation using 
a Fermi temperature of 5 meV. The s and p decomposition 
is shown for both types of atoms (B-I and B-II) in the ai2B 
structure. 



2.5 

3.3 
3.1 
2.9 
2-7 
2.5 
2.3 
2,1 
1.9 
1.7 
1.5 
1.3 
1 1 
0.9 
0,7 
0.5 
0.3 
0-1 
-0.1 
-0.3 




48.0 

46.0 

44,0 

42.0 

40.0 

38.0 

36.0 

34.0 

32,0 

30.0 ^ 

28.0 > 

26.0 A 

24.0 ^ 

22.0 O 

20.0 £ 

18,0 z 

16 IJJ 

14.0 

12.0 

10,0 

BO 

6.0 

40 

2.0 

00 
-2.0 
-4.0 



FIG. 3: The electronic band structure of fee B at the lattice 
constant 4.60 a.u. 



may have the same origin in both materials. 

The band structure of fee B confirms a metallic phase 
for all lattice parameters, with a rapidly increasing band 
width going to higher pressures. The occupied band with 
is 1.43 Ry at equilibrium and 1.88 Ry at a = 4.60 a.u. 
The energy bands for a = 4.60 a.u. are shown in Fig. |^. 
A fairly flat band appears near Ep in the XW direction 
at all lattice parameters. At X the band has strong s 
character, becoming p-like at W. This situation is remi- 
niscent of MgB2 where a similar flat band just above Ep 
has been identifijed,as the origin of superconductivity in 
this compound. crLratl This band appears to be responsi- 
ble for a peak in the DOS at Ep shown in Fig. ^. Also 
in Fig. we find the s, p, and d distribution of the DOS. 
In Fig.^ we confirm that the strongest DOS component 
is that from p-states. The values of DOS at Ep enter 
Eq. (|l|) above for the evaluation of the Hopfield parame- 
ter rj. The input and results using Eq. (fil) at four lattice 
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FIG. 4: The electronic density of states of fee B at the lattice 
constant 4.60 a.u., as determined by the tetrahedron method. 



parameters are listed in Table |. 

We note that in Table || the ^-components of the 
DOS are given within the MT spheres (as is the case 
in Fig. ^),and therefore they do not add up to the value 
of N{Ep). The p-component, Np{Ep), is, as expected, 
dominant and is about 60% of the value of the total in- 
side the MT. The p-d scattering, resulting from the £ ~ 1 
term in Eq. (Q), gives by far the largest contribution to 
the Hopfield parameter rj. For the high pressure case 
(a — 4.60 a.u.) rjpd = 0.79 x rjtot- The value of rjtot is 
approximately the S|ame as the value of 77 = 14.13 eV/A^ 
previously reportedu for metallic hydrogen at a higher 
pressure of 467 GPa. It is also interesting to note that the 
boron value of rj at equilib|rium is very close to the value 
77 = 7.627 eV/A^ for Nb,Ej a typical transition metal su- 
perconductor where the d-f scattering is the dominant 
term in Eq. (Q). 

To evaluateJiie critical temperature Tc we have used 
the McMillanEiJ approach defining an electron-phonon 
coupling constant A = rj/M <uj^>. Here 77 is calculated 
from Eq. (|l|) as discussed above. For the average phonon 
frequency <uj> we followed Eremets et ala and chose a 
range of values from 1200K to 1400K. The resulting A is 
in the range 0.53 to 0.39 respectively. Continuing we use 
the McMillan equation: 
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We solve Eq. (||) for five values of the Coulomb pseu- 
dopotential /i* — 0.09 to 0.13, and in the above range of uj 
and A values. The corresponding values of Tc are shown 
in Fig. ||. Within the uncertainty of the values of the 
frequency w, and near the value lo = 1250K our model 
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TABLE I: Radius of the muffin-tin sphere Rs, Fermi level _Ef, Total DOS at Ef, angular components of the DOS A'^, free- 
scatterer DOS N^^\ scattering phase shifts 5{£), and £ components of the Hopfield parameter ri and the total value of 77. The 
three columns are headed by the lattice parameters and the corresponding pressures. 
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FIG. 5: The superconducting transition temperature Tc of 
fee Boron, as determined by the McMillan equation (jij), as a 
function of the RMS phonon frequency <uj> at various values 
of fl*. 



predicts Tc values in the range of the experiment which 



validates the assertion of an electron-phonon mechanism. 

Finally, we want to comment on the pressure depen- 
dence of Tc- Our calculations show that the parameter 77 
goes from the value 7.82 eV/A^ at the equilibrium fee vol- 
ume to the value of 14.59 eV/A^ at a pressure of 307 GPa. 
It is important to note that N{Ep) follows the opposite 
trend, i.e., it reduces with increasing pressure. Assuming 
that the variation of <a;> with pressure is not as strong 
as the variation of 77 we propose that the rapid increase 
of 77 is responsible for the observed increase of Tc with 
pressure. 

In summary, we find that LAPW calculations, the 
RMT theory, and a McMillan analysis for Tc give a good 
description of superconductivity in boron at high pres- 
sures. 
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